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Abstract 

In supersymmetric models the lightest Higgs boson may decay with a sizable 
branching ratio into a pair of light neutralinos. We analyze such decays within the 
context of the minimal supersymmetric standard model with R-parity violation, where 
the neutralino itself is unstable and decays into Standard Model fermions. We show 
that the R-parity violating couplings induce novel Higgs decay signals that might fa- 
cilitate the discovery of the Higgs boson at colliders. At the LHC, the Higgs may 
be observed, for instance, through its decay -via two neutralinos- into final states 
containing missing energy and isolated charged leptons such as l^i^, £ ± £ ± , 3£, and 
41. Another promising possibility is the search for the displaced vertices associated 
with the neutralino decay. We also point out that Higgs searches at the LHC might 
additionally provide the first evidence of R-parity violation. 



1 Introduction 



The discovery of the Higgs boson is probably the most important goal of the LHC In the 
Standard Model, the Higgs branching ratios depend only on the unknown Higgs mass, which 
is constrained to be larger than 114 GeV. In the minimal supersymmetric Standard Model 
(MSSM) the Higgs sector is more involved, as it includes two Higgs doublets pQ. Yet, the 
mass of the lightest Higgs is very restricted, typically below 135 GeV. Another remarkable 
feature of the MSSM Higgs is that it may decay into two light neutralinos. Besides being 
invisible, such a decay causes a suppression of all other branching ratios of the Higgs boson, 
rendering more difficult its observation. In R-parity violating models, however, a new twist 
occurs. There, the neutralino is unstable and may decay within the detector into Standard 
Model fermions. Thus, the decay h — > \Xi with \ subsequently decaying into light particles, 
becomes visible and new signals for Higgs decays at colliders appear. That is exactly the 
situation we aim to study in this paper. 

The most general supersymmetric version of the Standard Model is phenomenologically 
inconsistent, for it includes lepton and baryon number violating operators that would induce 
fast proton decay. In the MSSM, an ad hoc discrete symmetry known as R-parity is imposed 
to prevent the decay of the proton. R-parity additionally implies the conservation of lepton 
and baryon number as well as the stability of the lightest supersymmetric particle -the LSP. 
Assuming R-parity, however, is not the only way of preventing proton decay. Lepton parity 
and baryon triality [2] are among the alternative discrete symmetries that forbid proton 
decay but allow for R-parity violation and, respectively, baryon or lepton number breaking 
couplings. R-parity, therefore, is not an essential ingredient of low energy supersymmetry. 

Supersymmetric models with broken R-parity are well motivated extensions of the Stan- 
dard Model and have been amply considered in the literature. They feature a rich phe- 
nomenology, markedly different from that of the MSSM. The LSP, for instance, is unstable 
and consequently it is no longer bound to be a neutralino; any supersymmetric particle 
can be the LSP [3]. And since the LSP decays into Standard Model fermions, distinctive 
decay patterns and collider signals are expected [U El [6] . R-parity violation might also be 
at the origin of neutrino masses and mixing [7J. The bilinear R-parity violating model, in 
particular, not only accounts for the observed values of neutrino masses and mixing angles 
[H [9] but it also predicts simple correlations between them and the LSP decay branching 
ratios |T0l [TTJ, IU El E] ■ These unique signals of R-parity violating models may soon be tested 
at the LHC as well as at future colliders. 

The mass of the lightest neutralino is not constrained by accelerator searches or precision 
experiments [T2| IT3] . If the GUT relation between gaugino masses is assumed, then the 
LEP limit on the chargino mass, m x ± > 100 GeV [TJ], translates into a lower bound on the 
neutralino mass: Mi ~ m x > 50 GeV[14|. If the assumption of gaugino mass unification is 
not made, however, there is no general limit on the mass of the lightest neutralino. Thus, 
neutralinos with masses below 50 GeV, light neutralinos, are perfectly compatible with 
present experiments. 

For the MSSM, the implications of such light neutralinos in Higgs boson decays were 
recently studied in [15J. There, after examining the dependence of the h — > XX branching 
ratio with the relevant supersymmetric parameters, it was shown that the decay into two 
neutralinos can even be the dominant decay mode of the Higgs boson, with a branching 
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Table 1: Strongest bounds on the R-parity violating couplings. 

ratio as large as 80%. Here, we extend such an analysis to supersymmetric models with 
R-parity violation. The R-parity violating couplings may cause the decay of the neutralino 
within the detector, transforming the former invisible decay h — > XX into a visible one. 

We will consider both bilinear and trilinear R-parity violating operators and study the 
different 3-body neutralino decays they induce. These decays give rise to novel Higgs decay 
signals that could be searched for at colliders. The most interesting searches are final states 
containing missing energy and isolated charged leptons such as £ ± £ T , £ ± £ ± , 3£, and 4£. In 
addition to these standard searches, the neutralino decay length can be long enough to leave 
a displaced vertex in the detector [H [TBI Ej- These searches may facilitate the discovery 
of the Higgs boson at the LHC and at the same time provide direct evidence of R-parity 
violation. 

The rest of this paper is organized as follows: In section [2] we discuss briefly the model 
considered as well as the constraints on the R-parity violating couplings. Section [3] is devoted 
to neutralino decays and presents our main results. The decay length is studied and several 
scenarios, in which different neutralino final states are present, are analyzed. Finally in 
section 0] we present our conclusions. 



2 The R-parity violating model 

The most general supersymmetric version of the Standard Model has a renormalizable su- 
perpotential given by 



W 



HmSSM + £ab 



k 5 



(1) 



where Wmssm is the lepton and baryon number conserving superpotential. In eq. (EE]), i,j, k 
run over the fermion generations, a, b are SU(2) indices, and a,/3,cr are color indices. The 
bilinear couplings i = 1,2,3, have mass dimension one and break lepton number. A 
and A' are dimensionless trilinear couplings that also break lepton number. There are 27 
independent \'^ k but only 9 \j k -they are antisymmetric in Baryon number is broken 
by the dimensionless couplings A^ fc , 9 of which are independent. In all, therefore, the R- 
parity violating superpotential contains 48 additional parameters. In addition there are the 
corresponding soft SUSY breaking terms. 
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To prevent proton decay, in the MSSM all the operators in ([I]) are forbidden by assuming 
the discrete symmetry known as R-parity. From a phenomenological point of view, however, 
the stability of the proton only requires that baryon and lepton number violating operators 
not be simultaneously allowed. And R-parity is not the only discrete symmetry able to 
ensure that. Baryon-triality and lepton-parity [2], for instance, are two well motivated 
symmetries that allow for either lepton or baryon number violating couplings, breaking 
R-parity but keeping the proton stable. 

Once R-parity is broken, the stability of the LSP is no longer guaranteed. The R-parity 
breaking terms in eq. (0Q) in fact induce the decay of the neutralino -mediated by a gauge 
boson or a scalar- to three Standard Model fermions. Such a decay constitutes the main 
difference between the MSSM and the R-parity violating model. 

Low energy data put a strong bound on some of the R-parity violating couplings in (pQ). 
For the bilinear couplings e^, the most stringent constraint comes from neutrino physics. The 
atmospheric mass scale can be generated at tree level if Y^i( e i v d + /ifi) 2 /Det(m x ) ~ m u /M 2 
where Det(m x ) is the determinant of the MSSM neutralino mass matrix, M 2 the SU(2) 
gaugino mass, Vi and Vd are the sneutrinos and H d vacuum expectations values. The solar 
mass scale, on the other hand, is induced only at the loop level and constrains the ratio 
|ej//i| < 10 -3 . For the trilinear couplings the strongest bounds come from double nucleon 
decay [IB] , neutron oscillations [l8j[T9], neutrino masses [20], [21], and neutrinoless double 
beta decay [22J. A complete and detailed list of the different constraints on the trilinear 
couplings can be found in [231 121] . The most important ones are summarized, for further 
use, in Table HJ Note here, that in principle one can rotate the "four vector" (Jl^ Lj) such 
that the bilinear terms are absent in the superpotential eq. ([I]) on the expense of changing 
the values of \jk and A- jfc (see e.g. [6]). The bounds in this table have to be understood in 
this particular basis. 



3 Neutralino decays 

The branching ratio h — > XX depends mainly on fi, tan/5, and m x [15], whereas the subse- 
quent neutralino decay is determined by m x , the R p violating couplings and gauge bosons 
or scalar masses. To study the decay of neutralinos originating in Higgs decays, we will 
consider a generic class of supersymmetric models featuring a non-negligible BR(/i — > xx) 
where all the parameters but the slepton and the squark masses are kept fixed. We take 

M x = 35 GeV, /i = 300 GeV, tan/3 = 5 , 
M 2 = M 3 = M A = 1 TeV , 

A = -1.7TeV, (2) 
nig > 800 GeV , 
m r > 200 GeV. 

The supersymmetric spectra thus obtained, illustrated in figure [3] satisfy the constraints 
from accelerator searches [H], the Higgs mass [25] . (g — 2) M [HI [26], and b — > |27j. To 
compute the spectrum and to evaluate the Higgs mass and other observables we use the 
FeynHiggs program [28]. According to it, lighter squarks -with all other parameters fixed- 
are not compatible with the LEP bound on the Higgs mass. 
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Figure 1: The typical supersymmetric spectrum we consider. 



Notice that for simplicity we assumed common soft-breaking masses for sleptons and 
squarks as well as a typical value of 1 TeV for M 2 , M 3 and Ma- Since Ma 3> Mz we are 
in fact working in the decoupling limit, where the interactions of the lightest Higgs boson 
become SM-like. 

Regarding tan f3, /x, and Ml, they were chosen so as to be in a region where BR(h — > xx) 
is non- negligible, and their values are rather typical within that region. Specifically, we get 
BR(h — > XX)— 21%, being bb the dominant decay mode with BR(/i — > bb)= 59%. As 
Mi <C M2,/i, the lightest neutralino is dominantly a bino, but it has a small higgsino 
component that generates the non-zero h\X coupling. 

If the neutralino decays outside the detector no new collider signals due to R-parity 
violating couplings are expected, as we would essentially be back to the MSSM case, where 
the decay h — > XX is invisible. We must, therefore, ensure that a significant number of 
neutralinos decays inside the detector. At the LHC, the Higgs is mainly produced through 
gluon fusion and its production cross section is huge -about 45 pb for our model. So, 
the neutralino decay length could be large and still yield a significant number of neutralino 
decays within the detector. This fact is illustrated in figure El where we quantify the fraction 
of neutralinos that decay inside a typical detector as a function of the decay length. This 
figure was generated with the PYTHIA program, version 6.414 [29], by taking into account 
only the production of the lightest CP-even Higgs. The MSSM parameters were varied 
according to eq. (jSJ) whereas a specific R-parity violating coupling was varied in the range 
[10~ 4 ,10 -1 ]. The points were selected by imposing the condition that the neutralino decays 
take place inside a cylinder of 3 m in the z direction and 0.9 m of radius, e.g. well inside the 
inner ATLAS or the inner CMS detector. From the figure we see, for instance, that if the 
proper neutralino lifetime is 10 m, about 7% of the decays will occur within the detector. 
Notice therefore that the predicted number of events will depend on the value of the R-parity 

1 The observed splitting between the two stop mass eigenstates is due to the non-zero value of the 
parameter A t 
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Figure 2: The relation between the fraction of neutralinos that decay within the detector 
and the neutralino decay length considering only neutralinos from Higgs decays. 

violating couplings. 

In what follows we will discuss the non-standard decays [13] of the Higgs boson that are 
induced by the R-parity violating couplings. 

3.1 Decays induced by bilinear terms 

Bilinear broken R-parity models are theoretically well motivated scenarios. They provide a 
simple framework that accounts for the observed values of neutrino parameters and that, in 
contrast with the seesaw mechanism, can be tested, through the decay properties of the LSP, 
in accelerator experiments. These models contain six lepton number violating parameters 
[30l [31] : €i, and their corresponding soft-breaking terms. These parameters are not entirely 
free, they are constrained by neutrino oscillation data. At present, the experimental data 
on neutrino oscillations indicates that [3"2"| I3"3"]: 

tan 2 12 = 0.47 ± 0.05, tan 2 6 23 = 0.83+°;^, sin 2 13 < 0.019 
Am 2 x = 7.67+g; 2 . 2 x 10" 5 eV 2 , Am 2 31 = 2.46 ± 0.15 x 10~ 3 eV 2 . (3) 

In our analysis, we will demand compatibility at the la level between these data and the 
six bilinear parameters. 

Neutralino decays in bilinear broken R-parity models are due to the mixing between 
neutralinos and neutrinos. The bilinear soft breaking terms, indeed, induce non-zero vevs 
for the sneutrinos that give rise to a mixing between leptons and gauginos and between 
sleptons and higgses. Such mixing allows the neutralino to decay -via a Z°, W^, sfermion 
or Higgs exchange- into the final states ViVjV^ Viqq, i^iltl^ or I qq'. 

Apart from these contributions induced by the mixing, there are additional ones due to 
the effective trilinear couplings |6j A233 = h T e 2 /n and A 333 = h b e 3 /fi. These new contribu- 
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Figure 3: The neutralino lifetime as a function of the slepton mass for different values of the 
bilinear R-parity violating couplings consistent with neutrino oscillation data at 1 a level. 

tions, mediated by slepton and squarks, give rise to neutralino decays into the final states 
utt, vfir, and vbb. 

Figure [3] shows the neutralino lifetime as a function of the slepton mass for models with 
bilinear R-parity violation. The supersymmetric spectrum was chosen according to equation 
(J2]) and the figure was generated with a private version of SPhencU El that includes bilinear 
R-parity violation. For any given value of the slepton mass, there is a spread in the neutralino 
lifetime that is due to the uncertainty on neutrino parameters. From the figure we see that 
the neutralino decay length has an upper bound of roughly 8 meters and that it is always 
larger than about 50cm. That means that, according to figure [21 between 7% and 40% of 
neutralino decays will occur within the detector. 

Two different regions can be easily distinguished from figure [3J For low slepton masses, 
m i ~ 800GeV, the neutralino lifetime increases with the slepton mass. Neutralino decays 
in this region are thus mediated by sleptons and induced by the effective trilinear couplings. 
For larger slepton masses, instead, the neutralino lifetime becomes essentially independent 
of rrig. In this region neutralino decays are mediated by gauge bosons and induced directly 
by the mixing. This picture is confirmed by figure HI where we show the corresponding 
neutralino branching ratios as a function of the slepton mass. As expected, the dominant 
decay modes at low slepton masses are r ± r T z/ and r ± / T z/ (/ = e, ji) whereas for large slepton 
masses several final states have sizeable branching ratios. 

2 This version can be obtained from W.P. 
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Figure 4: Neutralino decay branching ratios as a function of the slepton mass. The bilinear 
R-parity violating couplings are consistent with neutrino oscillation data at 1 a level. 



Models with bilinear R-parity violation have thus two remarkable features. On the one 
hand, the neutralino lifetime can be predicted with certain confidence. It lies between 50 
cm and about 8 meters; so it is rather large but it has a known upper bound. On the 
other hand, the neutralino decay products are rather uncertain and strongly depend on the 
sfermion masses. However, certain rations of branching ratios are predicted in terms of 
neutrino mixing angles In general, several final states with sizable branching ratios are 
expected. 



3.2 Decays induced by trilinear terms 

The R-parity violating couplings A, A', A" induce 3-body neutralino decays into Standard 
Model particles. To study such decays, we assume that, in turn, only one of the couplings, 
say A122, dominates and all others are negligible. First we consider the baryon number 
violating couplings A"- fe , and then the lepton number violating A^- fc and A^fc. For each case 
we compute the neutralino lifetime, function of the couplings and the sfermion 

masses. With that information, we determine the range of R-parity violating parameters 
that lead to neutralino decays inside the detector, find the new Higgs decay signals they 
induce, and briefly analyze the possibility of observing them at colliders. 

For simplicity, we work in the approximation where all final state fermions are massless^]. 
This approximation breaks down only if there is a top quark in the final state. Such a decay, 

3 We have checked that even taking properly into account changes the total width only slightly. 
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however, is not kinematically allowed. The neutralino decay width, then, has the generic 
form 



c,g 2 {\, A', A"} 2 



r = ' J -f(m?), (4) 

x 1536tt 3 Jy f h v > 

where Cf is the color factor, { } denotes one of the couplings, and fijnj) is a function of 
dimension m -4 that depends on the masses of the intermediate sfermions 



3.2.1 Decays induced by A" 

The trilinear coupling X"- k may induce the decay of neutralinos into 3-quark final states, such 
as uds and csb, leading to a Higgs boson that decays into a six-quark final state. Such decays 
were previously considered in [36J, albeit in a different scenario. Indeed, it was assumed 
in [36] that the Higgs boson had a mass around 100 GeV and had been missed by LEP 
searches because of its dominant decay into six quarks. Such a situation, however, is only 
possible in a very restricted portion of the viable parameter space. We, instead, consider a 
more generic framework where the Higgs is compatible with the usual LEP bound and has 
a non-dominant BR(/i — > xx)- 

Three different diagrams, respectively mediated by Ui, dj, and c4, contribute to the 
neutralino decay induced by a given \'[- k . The possible final states are Uidjd k (j ^ k) and 
its conjugate Uidjdk. Hence, the total decay width is simply given by T x = 2T(x Uidjdk). 
Besides A", T x will only depend on the squark masses. Figure [5] shows the neutralino lifetime 
as a function of the squark mass for A" = 10" 1 , 10~ 2 , 10 -3 , 10 -4 . The resulting neutralino 
lifetime varies, for trig < 1.6 TeV, approximately between 100 fim and 1 km. 

Notice that the decays induced by the couplings X'^ k are kinematically forbidden, as they 
give rise to final states containing a top quark. Moreover, due to the strong constraint that 
exist on A'/ 12 and A'/ 13 , see Table [fl neutralino decays into uds and udb are very suppressed 
and take place outside the detector. 

The unique signal from Higgs decays in this case is then a six-quark final state. But it is 
not known whether such signal could be observed over the QCD background. An interesting 
possibility, put forward in |37j, is the search -at LHCb- for the displaced vertices associated 
with the neutralino decay. As observed in figure [5j the neutralino decay length may be larger 
than 100/iin, leaving a displaced vertex; and the LHCb detector, in contrast to ATLAS and 
CMS, is well suited to study such events. 



3.2.2 Decays induced by A' 

Neutralino decays induced by the trilinear coupling A' contain leptons in the final state and 
are, therefore, easier to observe. An analysis related to ours was presented several years ago 
in [38] . They considered the special case of Higgs production through vector boson fusion 
and studied the Higgs signals induced by only certain trilinear couplings A and A' assuming 
gaugino mass unification. 

The coupling A- jfc gives rise to two different final states (plus their conjugates): 

X -> eiUjd k , (5) 
X -> Uidjdk. (6) 
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Figure 5: The neutralino lifetime as a function of the squark mass for different values of the 
R-parity violating couplings A". 

And each of them receives contributions from three different diagrams. The decay (jSJ) may 
have e-i, Uj, and dk as intermediate particles while Pi, dj and dk mediate the process (jSj). 
The total neutralino decay width is then given by 



In this case the neutralino decay width depends on A', the squark masses, and the slepton 
masses. We show, in figure El the neutralino decay length as a function of the slepton mass 
for A' = 10 _1 , 10~ 2 , 10 -3 , 10~ 4 and nig = 800 GeV. Notice that for large slepton masses, 
the diagrams mediated by squarks tend to dominate over those mediated by sleptons and 
consequently the curve becomes rather flat. From the figure we see that the neutralino 
lifetime varies, for < 1.6 TeV, between 0.1/xm and 100m. Given that the final state e{tdk 
is not kinematically allowed, the couplings X r i3k induce neutralino decays only into Vibdk- 
For those couplings, therefore, the neutralino lifetime is actually a factor of two larger than 
shown in figure 

The possible signatures of the decay of the Higgs boson due to the coupling A' are: 

1. Zero lepton, jets, and missing energy. 

2. One lepton, jets, and missing energy. 

3. Opposite sign lepton pair and jets. 

4. Same sign lepton pair and jets. 

Since standard searches for supersymmetry at colliders usually rely on missing energy 
signals, events with no missing energy, such as 2£+jets, might simply be rejected at the 




(7) 
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Figure 6: The neutralino lifetime as a function of the slepton mass for different values of 
the R-parity violating couplings A'. The common squark mass was set to 800 GeV. 

trigger level and never be available to study. They will, however, give rise to a displaced 
vertex provided that A' < 0.01. The decays induced by the couplings X' i3k always give rise 
to jets plus missing energy signals. Even after the degrading in missing energy compared 
with the case of stable neutralino, the signal with jets and neutrinos, which has at least a 
50% branching, has good potential to be discovered at the LHC [T7] . A generic expectation 
of this scenario is that the Higgs should be discovered at LHCb [37J. 

3.2.3 Decays induced by A 

The couplings A^ induce neutralino decays into final states containing two charged leptons 
and one neutrino. A given A^ might lead to two different final states (plus their conjugates): 

X -> CifjCk (8) 
X -> Pi&je k . (9) 

Hence, a neutrino is always present in the final state. As before, three different diagrams 
contribute to each final state and the resulting decay width simply scales as l/m~. Figure 
[7] shows the neutralino decay length as a function of the common slepton mass for different 
values of the coupling A. 

The two signatures with missing energy we mentioned in 13.2.21 -numerals 1 and 2- will 
also be present in this case, though the jets come from the hadronic decay of the tau lepton 
and not from final state quarks. Additionally, the A couplings also give rise to new signatures 
with two or more leptons. They are 

1. Opposite sign lepton pair, jets, and missing energy. 
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Figure 7: The neutralino lifetime as a function of the slepton mass for different values of 
the R-parity violating couplings A. 

2. Same sign lepton pair, jets, and missing energy. 

3. 3 leptons + jets + missing energy. 

4. 4 leptons + missing energy. 

Notice that the decays induced by A, in contrast with those due to A', always lead to missing 
energy -from the final state neutrino. These new signatures with two or more leptons are 
particularly significant because thanks to their low backgrounds they have good chances to 
be discovered at LHC [THEE]. 

4 Conclusions 

We studied the decay of the Higgs boson into Standard Model particles within the context 
of the MSSM with R-parity violation. The decay proceeds via h — ► XX followed by the R- 
parity violating neutralino decay into light fermions. We pointed out that neutralino decays 
induced by the trilinear R-parity violating couplings may occur inside the detector and, as 
a result, give rise to novel Higgs decay signatures that may facilitate the discovery of the 
Higgs boson at the LHC. Particularly appealing -because of their low backgrounds- are the 
decays into final states containing three or four leptons and missing energy. Such decays 
are caused by the lepton-number violating couplings A^ and could be easily observed at 
the LHC. Another promising possibility is the observation of displaced vertices associated 
with the neutralino decay. Thus, the discovery of the Higgs boson at the LHC might also 
provide direct evidence of R-parity violation. 
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